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Abstract 



Negative ion detachment in bichromatic laser field is considered within 
the adiabatic theory. The latter represents a recent modification of the fa- 
mous Keldysh model for multiphoton ionization [L. V. Keldysh, Zh. Eksp. 
Teor. Fiz. 47, 1945 (1964) [Sov. Phys.-JETP 20, 1307 (1965)]] which makes 
it quantitatively reliable. We calculate angular differential detachment rates, 
partial rates for particular ATD (Above Threshold Detachment) channels and 
total detachment rates for H~ ion in a bichromatic field with 1 : 3 frequency 
ratio and various phase differences. Reliability of the present, extremely sim- 
ple approach is testified by comparison with much more elaborate earlier 
calculations. 

PACS numbers: 32.80.Rm, 32.80.Fb 
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I. INTRODUCTION 



Interest to the photoionization of atoms in a trichromatic laser field both in theory (see, 
for instance, Ref. [I | [ 13 1 ) and in experiment |l j [11] [16|] seems to stem first of all from 
the effect of the phase control, i.e. dependence of the observables on the difference of field 
phases <p. 

The calculations have been carried out previously for ionization of the hydrogen atom 
in two laser fields with a frequency ratio 1 : 2 M , 1 : 3 0] and 2:30. Potvliege and Smith 
H presented results for various frequency ratios and initial states. Different schemes have 
been employed, but all of them imply numerically intensive work. 

For the multiphoton electron detachment from negative ions some analytical treatment 
exists || which aims to investigate qualitative features of the process, mostly in the 
case when one or both fields are weak. The presence of large number of parameters in the 
problem sometimes makes results of analytical studies not directly transparent. Quantitative 
reliability of these approaches has never been assessed. This situation looks particularly 
unsatisfactory since the multiphoton electron detachment from negative ions presents unique 
situation when quantitative results can be obtained by analytical method in a broad range 
of parameters characteristic to the problem. Indeed, it has been demonstrated recently 
by Gribakin and Kuchiev (T7| that proper application of the well-known Keldysh (TP] 
model to multiphoton detachment provides an extremely simple scheme that gives very 
reliable results for the total rates as well as for ATD (Above Threshold Detachment) spectrum 
and ATD angular distributions. This adiabatic approximation ensures an accuracy which is 
comparable with that of the most elaborate numerical developments and works unexpectedly 
well even outside its formal applicability range, i.e. even for small number n of photons 
absorbed. The evidences of good performance of the Keldysh model for the total rates were 
presented also in the earlier paper [^5[ . 



Recently the adiabatic approach was extended by the present authors [^6J to the case 
of bichromatic field. The practical applications were carried out for the case of frequency 
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ratio 1 : 2 when in addition to the phase effects another unusual phenomenon exists, namely 
the polar asymmetry of the angular distribution of photoelectrons. Unfortunately no other 
quantitative data for photo detachment in this case is available which makes comparison 
impossible. 

The main objective of the present study is to assess quantitatively an accuracy of the adi- 
abatic scheme by comparison with the previous calculations carried out by Telnov et al Jll| 
in case of 1 : 3 frequency ratio. For this ratio the polar asymmetry is absent, but the phase 
effects persist. The calculations by Telnov et al JTTJ are based on sufficiently sophisticated 



numerical scheme providing a useful benchmark. We present (Sec. |J) complete comparison 
of the results by considering angular differential detachment rates, heights of ATD peaks 
and total detachment rates. It should be emphasized that the angular differential rates are 
most sensitive to the formulation of the model representing an ultimate test for the theory, 
as discussed in Sec.|T|. We draw also some general conclusion on the relation between the 
adiabatic approach and the numerical calculations within the one-electron approximation. 

II. RESULTS 



The adiabatic theory of two-color detachment was outlined in our previous paper [26 
where the reader can find all the details of calculation. Here we only write down the expres- 
sion for the electric field strength in the bichromatic laser field with 1 : 3 frequency ratio in 
order to specify the definition of the field phase difference <p 

F{t) = Fx cos cut + F 2 cos(3ujt + <p) . (1) 

Fi, F 2 are the amplitude vectors for the fundamental frequency u and its third harmonics 



respectively. Below we consider, just as in Ref. [11|, the case when both fundamental field and 



its third harmonics are linear polarized with Fi \\ F 2 . Then the differential photoionization 
rate depends only on the single angle 9 between the photoelectron translational momentum 
p and the vectors F\ \\ F 2 . Atomic units are used throughout the paper unless stated 
otherwise. 



Our calculations for H~ detachment are carried out for the parameters of H~ as before 
(k = 0.2354, A = 0.75). We choose two sets of field intensities l\ and I2 for the 
fundamental frequency to = 0.0043 (that of CO2 laser) and its third harmonics, same as 
in the paper by Telnov et al |IT| , namely (i) I\ = 10 10 W/cm 2 , I2 = 10 9 W/cm 2 and (ii) 
I x = 10 10 W/cm 2 , J 2 = 10 8 W/cm 2 . 

In case of the frequencies ratio 1 : 3 considered here the field ([!]) does not possess polar 
asymmetry (i.e. asymmetry under inversion of the z axis directed along F\ || F2). Therefore 
the differential detachment rate does not change under the transformation 9 71 — 9. This 
allows us to show plots only for > 9 > domain. 

Figs. [l|-|6| show the differential detachment rate as a function of the angle 9 for three 
lowest (open) ATD channels and for two sets of field intensities. The system Hamiltonian is 
a 27r-periodic function of the phase parameter ip. We show our results for ip = 0, ip — ±| and 
(f = tt. The transformation ip =$> n — <p leaves the Hamiltonian invariant only if t is replaced 
by —t. As stressed in Refs. [Q], the problem is invariant under the time inversion operation 
provided the final-state electron interaction with the atomic core is neglected. This is the 
case in the present model. Therefore our differential ionization rates are the same for ip and 



—ip. The calculations by Telnov et al [0 do take into account the final state electron-core 
interaction. Therefore they show some difference between the angular differential rates for ip 
and —<p. However, it proves to be quite small for low ATD channels as seen from the plots. 
The importance of the interaction between the emitted electron and the core has been 



first pointed out by one of the present authors [p^] . In this paper several phenomena has been 



predicted for which this interaction plays crucial role. The related mechanism was named 
"atomic antenna". In the recent literature the final state interaction is usually referred to 
as rescattering. In our problem the rescattering effects are enhanced for high ATD channels 
as discussed below. 

The results of our extremely simple theory are compared in Figs. [1|-|6] with the previous 
numerical calculations by Telnov et al |TT]] which are rather involving. Being carried out in 
the one-electron approximation, they employ an accurate model for the effective one-electron 



potential in H |28| , complex-scaling generalized pseudospectral technique |29| to discretize 



and facilitate the solution of the time-independent non-Hermitian Floquet Hamiltonian for 
complex quasienergies and eigenfunctions, and calculation of the electron energy and angular 



distributions by the reverse complex-scaling method j30] . As a lucid illustration of simplicity 
of the present approach it is worthwhile to stress that it does not rely on any particular form 
of an effective one-electron potential albeit employs only two parameters k and A governing 
the asymptotic behavior of the initial bound state wave function. 

From Figs. [l|-|6] one can see that the adiabatic approximation ensures good quantitative 



agreement with calculations by Telnov et al [11|. In particular, positions of maxima and 



minima in the angular photoelectron distribution are well reproduced. This demonstrates 
that the adiabatic approach correctly describes the nature of the structure as due to inter- 
ference between the electron waves emitted at various (complex-valued) moments of time. 
Indeed, within the adiabatic theory 0] the ionization amplitude is expressed as a sum of a 
number of interfering contributions. Mathematically they come from different saddle points 
in the approximate evaluation of the integral over time that emerges in the Keldysh (TP| 



model. Physically they correspond to the coherent emission of photoelectron at different 
moments of time. For our particular frequency ratio 1 : 3 the sum contains 6 interfering 
contributions as compared with 4 term for 1 : 2 frequency ratio [2(| and 2 terms for one- 



color detachment [[L7j] ||18|| . Generally this suggests that in the former case more complicated 
angular patterns emerge. Probably one can find here a correlation with an alternative in- 
terpretation in the multiphoton absorption framework. The latter argues [[Ll|] that the 
angular distribution structure in 1 : 3 case is more complicated than for 1 : 2 ratio since all 
the pathways leading to a continuum state with the same energy interfere in the 1 : 3 case 
whereas a considerable pattern of non-interfering pathways exists for the the 1 : 2 case due 
to parity or energy restrictions (each pathway is characterized by the number of photons of 
different colors absorbed successively). 

The partial detachment rate for each ATD channel are shown in tables | and [TJ for two sets 
of field intensities. The agreement is good for low ATD channels; note that the rescattering 



effects which generate dependence on the sign of tp are manifested in the partial rates even 
less than in the angular distributions shown in Figs. [1]-|6]. For higher ATD channels with low 



rates the difference between the present results and those of Telnov et al |TT| becomes more 
pronounced. This behavior could be interpreted as increasing importance of rescattering for 
high ATD peaks. The manifestations of this effect were observed recently in experiment ET 
and are currently vividly discussed in the literature |27j [^] [|33|]. 



III. CONCLUSION 

As a summary, the adiabatic approach provides quantitatively reliable tool for investi- 
gating two-color photodetachment of negative ions. In particular, the interference structure 
in the photoelectron angular distributions as well as the phase effects are correctly described. 
Since generally the interference phenomena are known to be most sensitive to the details 
of theoretical description, one can conclude that the present theory had successfully passed 
the stringent test. 

The Keldysh scheme |19| is known to be gauge-noninvariant. Importantly, the calcula- 



tions within the adiabatic approach |l7j [^6j employ the dipole-length gauge for the laser 



field thus stressing contribution of the long-range asymptote of the initial bound state wave 
function. The use of the length gauge together with the adiabatic approach (i.e. integration 
over time by the saddle point method, see Refs. fl7j [0 |2(J and discussion in Sec.|J) ren- 
der self-consistent character to the theoretical scheme. Indeed, the exact evaluation of the 
integrals does not add to the accuracy of the result as compared with the use of the saddle 
point method. This is because in the former case the integral absorbs the contributions from 
the wave function outside its asymptotic domain, where in fact it is known with much lower 
accuracy (being, in particular, influenced by the effects beyond the single active electron 
approximation) . 

The method is straightforwardly applicable to the negative ions with the outer electron 
having non-zero orbital momentum, such as halogen ions, which could be easier accessible 
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for the experimental studies (for the one-color detachment such applications could be found 
in Ref. ||17|| ). Technically the calculations within the adiabatic approach are extremely 
simple reducing to finding the roots of polynomial and substituting them into an analytical 
expression |26| (the related Mathematical 0] program takes only few lines). It should be 
recognized that the single active electron approximation itself introduces some intrinsic error. 
It seems that often this error could be comparable with the difference between the result of 
numerical one-electron calculations and these of the adiabatic approximation. Uncertainty 
of the one-electron approach in principle could be removed within the two-electron approach 
which however consumes much more efforts. The two-electron calculations which has been 
carried out recently show that the one-electron approximation is generally sufficient unless 
one is particularly interested in the subtle resonance effects |25| |35| (the calculations 
beyond one-electron approximation are currently possible only for small number of absorbed 
photons). For high ATD channels with low intensities the adiabatic approximation becomes 
less reliable due to increasing role of rescattering effects neglected in the present form of 
the approximation. It seems however that relatively simple modifications of the adiabatic 
approximation could be carried to include rescattering effects. 

Reliability of the results obtained above for the simple one-electron problem with rescat- 
tering neglected is highly important in perspective, since they are to be included as a con- 
stituent part in the treatment of much more sophisticated one-electron and many-electron 
problems governed by the antenna mechanism [p7[] p7| p8| . 
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TABLES 

TABLE I. Partial rates for the H~ detachment by the laser wave with the frequency to = 0.0043 
and its third harmonics with the intensities respectively l\ = 10 10 W/cm 2 and I2 = 10 9 W/cm 2 . 
The number of absorbed photons n refers to the fundamental frequency. In each block the upper 



figure gives present result and the lower one the result obtained by Telnov et al [11]. The number 
in square brackets indicate the power of 10. 



One-colour Two-colour One-colour 

n fundamental harmonic 

ip = cp = 7T ip = \ir ip = —\k 



0.67[-9] 
0.72 [-9] 



0.47[ 
0.421 



7] 
7] 



0.54[- 
0.58[- 



0.22 
0.20 



-7] 

:-7] 



0.22 
0.21 



7] 

-7] 



0.20[-9] 
0.20[-9] 



0.11 
0.101 



-7] 
-7] 



0.23[ 
0.23T 



0.80[ 
0.711 



-8] 
-8] 



0.80[ 
0.731 



-8] 
-8] 



0.46[-10] 
0.30[-10] 



10 



0.41 [-10] 
0.39[-10] 



0.27[ 
0.261 



8] 
8] 



0.34[~ 
0.27[- 



0.39 
0.35 



0.39 

0.30 



8] 
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0.50[-ll] 
0.40[-ll] 



0.65[ 
0.721 



9] 
9] 



0.23[ 
0.16[- 



0.17[ 
0.151 



-8] 
-8] 



0.17[ 
0.101 



-8] 
-8] 



12 



0.74[-12] 
0.71 [-12] 



0.16[ 
0.20[ 



9] 
9] 



0.10[ 
0.71[ 



0.62| 
0.581 



9] 
9] 



0.62| 
0.30[ 



9] 
9] 



0.66[-13] 
0.86[-13] 



12 



0.21 
0.33 

0.64 
0.14 

0.16 
0.47 

0.36 
0.14 

0.68 
0.35 



-12] 
-12] 

-13] 
-12] 

-13] 
-13] 

-14] 
-13] 

-15] 
-14] 



0.38 
0.53 

0.88 
0.14 

0.20 
0.32 

0.52 
0.72 

0.15 
0.22 



-10] 
10] 



0.36[-9] 
0.27[-9] 



0.20[-9] 
0.20 [-9] 



0.20[-9] 
0.85[-10] 



-11] 0.11 [-9] 0.58 

-10] 0.97[-10] 0.69 

-11] 0.32[-10] 0.16 

-11] 0.33[-10] 0.23 

-12] 0.82[-ll] 0.44 

-12] 0.12[-10] 0.78 



-12] 
-12] 



0.21[-11] 
0.43[-ll] 



0.12 
0.27 



-10] 
-10] 

-10] 
-10] 

■11] 
-11] 

-11] 
11] 



0.58 
0.32 

0.16 
0.17 

0.44 
0.10 

0.12 
0.60 



-10] 
-10] 

-10] 
-10] 

-11] 

-10] 

-12] 
-11] 



0.16[-15] 
0.31[-15] 



Total 



0.92[-9] 
0.96[-9] 



0.62 [-7] 
0.56[~7] 



0.15[-7] 
0.14[-7] 



0.36 [-7] 
0.33[-7] 



0.36[-7] 
0.33[-7] 



0.46 [-7] 
0.30[-7] 
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TABLE II. Same as in table 1, but for the intensities h = 10 10 W/cm 2 and I 2 = 10 8 W/cm 2 . 



One-colour Two-colour One-colour 

n fundamental harmonic 

<f = If = 7T if = ^7T (f = — ^VT 



0.67[-9] 
0.72 [-9] 



0.54[-8] 
0.53[-8] 



0.41 [-9] 
0.36[-9] 



0.25[-8] 
0.24 [-8] 



0.25[-8] 
0.25 [-8] 



0.20 [-9] 
0.20 [-9] 



0.92[-9] 
0.93[-9] 



0.16[-9] 
0.18[-9] 



0.69 [-9] 
0.68[-9] 



0.69 [-9] 
0.68[-9] 



0.46[-13] 
0.30[-13] 



0.41[ 
0.39[ 

0.50[ 
0.40[ 

0.74[ 
0.71[ 

0.21[ 
0.33[ 

0.64[ 
0.14f 



-10] 
-10] 

-11] 
-11] 

-12] 
-12] 

-12] 
-12] 

-13] 
-12] 



0.23[-9] 
0.25[-9] 



0.54[ 
0.59| 

0.11! 
0.14[ 

0.21| 
0.45[ 

0.39[ 
0.211 



-10] 
-10] 

-10] 
-10] 

-11] 
-11] 

-12] 
-11] 



0.82 
0.82 

0.74 
0.65 

0.38 
0.33 

0.13 
0.13 

0.39 
0.47 



-io] 
-io] 

-io] 
-io] 

-io] 
-io] 

-io] 
-io] 

-11] 
-11] 



0.25 [-9] 
0.25 [-9] 



0.25 [-9] 
0.22 [-9] 



0.88[ 
0.88[ 

0.28[ 
0.28[ 

0.79[ 
0.82[ 

0.20[ 
0.241 



-10] 
-10] 

-10] 
-10] 

-11] 
-11] 

-11] 
-11] 



0.88[ 
0.66[ 

0.28[ 
0.22[ 

0.79[ 
0.94[ 

0.20[ 
0.49f 



-10] 
-10] 

-10] 
-10] 

-11] 
-11] 

-11] 
-11] 



0.66[-17] 
0.98[-17] 
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15 0.16[-13] 0.66[-13] 0.96[-12] 0.46[-12] 0.46[-12] 0.16[-21] 
0.47[-13] 0.11[-11] 0.17[-11] 0.70[-12] 0.25[-ll] 

16 0.36[-14] 0.11[-13] 0.21[-12] 0.99[-13] 0.99[-13] 
0.14[-13] 0.52[-12] 0.60[-12] 0.21[-12] 0.12[-11] 

17 0.68[-15] 0.17[-14] 0.42[-13] 0.21[-13] 0.21[-13] 
0.35[-14] 0.22[-12] 0.22[-12] 0.64[-13] 0.50[-12] 

Total 0.92[-9] 0.66[-8] 0.79 [-9] 0.36[-8] 0.36[-8] 0.46[-13] 

0.96[-9] 0.66[-8] 0.74[-9] 0.35[-8] 0.35[-8] 0.30[-13] 
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FIGURES 

FIG. 1. Detachment of H~ ion in trichromatic field with the frequencies uj = 0.0043 and 3uj 
and intensities I\ = 10 10 W/cm 2 and I2 = 10 9 W/cm 2 respectively. Differential detachment rate 
(in units 10 _8 a.u.) as a function of the electron emission angle 9 is shown for the first ATD peak 
(corresponding to absorption of n = 8 photons of frequency uj) and various values of the field phase 
difference tp as indicated in the plots. Open symbols show results of calculations by Telnov et al 
fll"T|| (in the tp = ±|vr plot the open circles show the results for ip = ^ir and open triangles these for 
(p = — ^7r). Solid curves show results of the present adiabatic theory (which coincide for ip = ^ir 
and ip = — \tx as discussed in the text). 

FIG. 2. Same as in Fig.l, but for the second ATD peak (corresponding to absorption of n = 9 
photons of frequency u). 

FIG. 3. Same as in Fig.l, but for the third ATD peak (corresponding to absorption of n = 10 
photons of frequency uj). 

FIG. 4. Same as in Fig.l, but for intensities I\ = 10 10 W/cm 2 and I2 = 10 8 W/cm 2 ; the 
differential detachment rate is shown for the first ATD peak (corresponding to absorption of n = 8 
photons of frequency uj). 

FIG. 5. Same as in Fig. 4, but for the second ATD peak (corresponding to absorption of n = 9 
photons of frequency uj). 

FIG. 6. Same as in Fig. 4, but for the third ATD peak (corresponding to absorption of n = 10 
photons of frequency uj). 
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